Dynamics of one-dimensional Bose liquids: 
Andreev-like reflection at Y-junctions and absence of the Aharonov-Bohm effect 
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We study one dimensional Bose liquids of interacting ultracold atoms in the Y-shaped potential 
when each branch is filled with atoms. We find that the excitation packet incident on a single 
Y-junction should experience a negative density reflection analogous to the Andreev reflection at 
normal-superconductor interfaces, although the present system does not contain fermions. In a 
ring interferometer type configuration, we find that the transport is completely insensitive to the 
(effective) flux contained in the ring, in contrast to the Aharonov-Bohm effect of a single particle in 
the same geometry. 

PACS numbers: 03.75.Kk 



Recently, guiding of atoms in low-dimensionally mag- 
netic traps has been actively studied 0]. Such systems 
provide an ideal opportunity to study coherent quantum 
dynamics of interacting many-body systems. The the- 
ory of one-dimcnsionally trapped atoms in equilibrium 
0, H, EL HI > as well as questions related to coherent dy- 
namics and nonequilibrium phenomena [|| , has been vig- 
orously studied. In this Letter we analyze a simple yet 
nontrivial example of the real-time dynamics of Bose- 
Einstein condensates (BEC) in the Y-shaped potentials 
(beam-splitters) Q, and in the related "ring interferom- 
eter" type geometry [§]. 

In contrast to previous works where the dynamics of 
wave packets in otherwise empty traps is studied, here 
we consider the situation where each one-dimensional 
branch of the Y-shaped potential is filled with interacting 
bosonic atoms, i.e. each branch contains one-dimensional 
Bose liquid. The Tomonaga-Luttinger (TL) liquid the- 
ory [!, H[ is a powerful method to study this problem. 
The whole system in the Y-shaped potential may be re- 
garded as three TL liquids connected at a junction, as in 
Fig. [TJ While there are parallels with the corresponding 
electron problem nj, as we will show in the follow- 
ing, a quite different generic behavior is expected in the 
Bose liquid. In the present case, instead of the conduc- 
tance, the real-time dynamics could be directly observed 
in experiments. Furthermore, as an application of the 
above analysis of the single Y-junction problem, we also 
consider the transport in a ring type interferometer ge- 
ometry. 

We have two main results in this Letter. The first 
is the prediction of the negative density reflection at a 
junction (Fig. [1]). In the discussion of Y junctions for 
systems of interacting electrons, such negative reflection 
has been interpreted as the Andreev reflection, which oc- 
curs in metal-superconductor interfaces due to Cooper 
pair formation. On the other hand, the negative reflec- 
tion at the BEC Y-junction can be understood rather 



simply using an analogy with electromagnetic transmis- 
sion line three-way junction. For an ideal three-way junc- 
tion, there should be no voltage difference between the 
connected ends of the three lines, and there is conserva- 
tion of the total charge in a pulse. These two conditions 
imply that when the amplitude of an incident signal is 
one, signals transmitted into the other two lines will have 
amplitudes of 2/3 and the reflected signal will have the 
amplitude of —1/3. As we discuss below the same phe- 
nomenon should appear in a BEC Y-junction. 

The other main result is the absence of the Aharonov- 
Bohm (AB) effect in the ring type interferometer, which 
consists of two Y-shaped junctions as in Fig. [2] As a 
typical example of the AB effect, when the ring is empty, 
a single particle injected from one lead to the ring does 
not transmit to the other lead when the (effective) "mag- 
netic flux" inside the ring is the half unit flux quantum. 
However, we find that the transport between the leads in 
our system is completely insensitive to the magnetic flux, 
in the low energy limit. These two findings are manifes- 
tation of the collective nature of the Bose liquid, rather 
than single-particle physics. 

We start from the following one-dimensional Hamilto- 
nian 12j with the low-energy scattering characterized by 



the low-energy s-wave scattering: 



dx 



U 



2m 



d x ipt(x)d x tpj{x) + —Pj{x) 



(1) 



where we take an unit of h = 1. vjj(x), Pj{x) = 
tpj(x)x/)j(x) are, respectively, the annihilation and the 
density operator of the boson with the mass m, on the 
j-th branch, (j = 1,2,3.) The effective interaction 
strength U > is determined by the s-wave scattering 
length. We take x > 0, so the physics at the junction is 
represented by the boundary condition (b.c.) at x = 0. 
In case of the junction problem, the b.c. is a non-trivial 
problem and plays a crucial role. First of all, we have 
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to identify what kind of b.c. describes the physics in the 
low-energy limit. 

Any b.c. should satisfy the total current conservation 

Jx{t,x= 0) + J 2 (t,x = 0) + J 3 (t,x= 0) = (2) 

where Jj(x) is the current operator of the Bose liquid in 
the j-th branch. The simplest possible b.c. corresponds 
to an infinitely strong barrier at the junction, which pro- 
hibits any particle to be transmitted through the junc- 
tion. It would give Jj(t,x = 0) = for j = 1,2,3. 
Naturally, this satisfies the current conservation ([2]). 

In practice, there would be a finite tunneling ampli- 
tude r between the different branches. The effect of the 
tunneling can be included by adding the perturbation 



Ti-B = — 



= 0)Vi-i(a = 0) + h.c. 



(3) 



at the boundary. We are interested in the low-energy 
behavior of the junction in the presence of the tunneling. 

In order to study the low-energy physics, it is conve- 
nient to use the "bosonization" technique [9]. In terms 
of the TL boson field 9j(x) and its dual ifij{x), the Bose 
atom annihilation operator and density operator on the 
j-th branch are, respectively, represented as tpj(x) ~ 

[pC°> + BxOjix)/*] 1/2 e-^<*), Pj {x) ~ + ft^CaOAr, 
where we have retained only the leading terms, p^ rep- 
resents the expectation value of the atom density in the 
ground state. The low-energy effective Hamiltonian for 
equation (fT]) is given by 



Ho = 



3 

E* 



f dx^d^f+g-'idj,) 2 ] 
/7r Jo 



(4) 



The fields 9, (p satisfy the commutation relation 
[0{x), y>{y)] — ittH(x — y), where H(x — y) is the Heav- 
iside step function, v is the velocity of the collective 
modes, and g is the so-called TL parameter, in which the 
interaction is essentially contained. 

The TL parameter, for the interacting fermions i.e. 
electrons, usually satisfies g < 1 due to the Coulomb re- 
pulsion. However, g > 1 occurs rather naturally in a Bose 
liquid with a short-range repulsive interaction. For exam- 
ple, for the Bose liquid with a (^-function interaction [l3j |. 
the TL parameter g can be determined from the Bethe 
Ansatz exact solution. It turns out that g is always larger 
than 1 for this model, interpolating between g = oo in 
the non-interacting limit and g = 1 in the strongly inter- 
acting (or dilute) "Tonks gas" limit. [1,0, @| While g can 
also be less than 1 with other types of the interaction 14 1 , 
g > 1 would be rather generic for the interacting Bose 
liquid. 

In terms of the bosonization, the "disconnected" limit 
ipj(0) = corresponds to the Dirichlet boundary condi- 
tion dt0j(0) = for j = 1,2,3. It is equivalent to the 



Neumann b.c. on ipj. The tunneling perturbation ([3]) is 
bosonized as H B ~ Fe~ l(l ^ (0) ^- l(0)) , where only the 
most relevant term is retained. The scaling dimension of 
this perturbation is given by g , and it is relevant when 
g > 1. Namely, the effective tunneling amplitude T grows 
as the energy scale such as the temperature is lowered. 
It is expected that F eventually grows to infinity in the 
low-energy limit, when g > 1. [151 ] 

In the strong coupling limit, the tp fields should be 
pinned as t/?i(i,0) = (^2(^,0) = f^{t, 0), because the tun- 
neling term acts as an infinitely strong potential on (p 
at the junction. This also implies dtfi = dty>2 = dtfs- 
Translating it into the b.c.s on 9j, we find 



J2d t 6 j (t,x = 0)=0, (5) 

d x d i (t,x = 0) = d x e j {t,x = 0) {i,j = 1,2,3), (6) 

where the first line comes from the current conserva- 
tion ((2]). The second line indicates that the density is 
continuous at the junction, reflecting that the junction 
is at the strong coupling limit. Using the boundary con- 
formal field theory techniques, we can see that all the 
possible perturbations to the above b.c. are irrelevant if 
g > 3/4 [13, EH in the case of Bose liquid. This guar- 
antees the stability of the "strong coupling" fixed point, 
supporting that T grows to infinity in the low-energy limit 
when g > 1. 

Now we discuss the dynamics of the system with the 
stable, strongly coupled b.c. Let us suppose that the 
density is altered locally from the equilibrium. The vari- 
ation of the density propagates in the system, and even- 
tually would be transmitted to other branches. 

We start from the Heisenberg equation of motion for 
the TL boson field 9j and dual field <pj (j — 1,2,3): 
d t 6j(t,x) = -vgd x ip {t,x), d t ifj{t,x) = -^d x 6j(t,x). 
Let us focus on the expectation values of the density and 
current. Instead of operators 9j and (fj , we then just need 
to determine the expectation values (9j) and ((fj), which 
obey the simple one-dimensional classical wave equation. 
The initial conditions are given as the expectation values 
of current and density at t = 0. 

For simplicity, we consider the situation where a left- 
moving packet is injected from one side (branch 1). Since 
the left movers go towards the boundary (x = 0) , and the 
right movers go away from the boundary, the initial con- 
ditions may be given as {pf(t = 0, x)) — +T>o(x)5ij , 
(pf(t = 0,x)) = p {0 \ where V = d x {9[{t = 0,x))/ir is 
the scalar which gives the initial condition on the expec- 
tation value of the densities. Solving the classical wave 
equation with the initial conditions, we obtain: 

(pi(t,x)) = pW+Vo(v + )-lvo(v-)H( V -), (7) 

(P2, 3 (*,z)) = p^ + ^Vo(rr)H(rr), (8) 
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FIG. 1: Reflection and the transmission of the injected cur- 
rent. The left panel describes the initial situation, and the 
right panel shows the situation after the reflection and the 
transmission. 



where if 1 = vt ± x. The second term in the rhs of 
the equation ([7]) represents the incident packet density, 
and the third term implies the (negative) reflection at 
the boundary. Namely, if a bump in the density moves 
through the branch 1 to the junction, a dip in the density 
is reflected back, (see FigQ) 

Eq. (|8|) shows the transmission density from the branch 
1 . Defining the transmission tensor as the ratio of the in- 
cident current and transmission currents, the transmis- 
sion tensor from the branch 1 to the branch 2 (branch 3), 
22,i (73,i), are obtained as T2,i = T31 = |. This means 
that the sum of the densities transmitted to branches 2 
and 3 is greater than the incident density from the branch 
1. The current conservation law is satisfied with the "neg- 
ative reflection" current in the branch 1, which appears 
as the second term in equation ([7]). The time evolution 
of the density ([7])-([8]), including the remarkable negative 
reflection, could be observed in experiments, with present 
techniques on atomic BEC. 

A similar behavior was found in the Y-junction of one- 
dimensional electron systems, in terms of the conduc- 
tance. It was argued as a consequence of Andreev re- 
flection similar to that in a superconductor-normal in- 
terface. What is surprising here is that the Andreev-like 
negative reflection is also found for the boson systems. 
The system of bosons is usually not related to Cooper 
pair formation, which is the underlying mechanism of the 
Andreev reflection. These results may be also derived 
from a linearized mean-field theory which is applicable 
to the weakly interacting regime. While the mean-field 
approximation breaks down in the strongly interacting 
regime, our derivation based on the TL liquid theory are 
valid even in the presence of strong fluctuation. Thus our 
results represent universal feature of general repulsively 
interacting boson systems. 

Let us extend our analysis to a "ring interferometer" 
type configuration, which consists of two Y-junctions. 
(see Fig E} We take the TL fields on the upper and lower 
arc in the ring to 6 u j , respectively. Then the symmetric 
combination 9 S = u +8i is coupled with the right and left 
leads, while anti-symmetric combination 6 a — 6 U — 81 is 
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FIG. 2: The schematic figure of the ring interferometer made 
of two Y junctions, with left and right leads, where 4> m an 
effective "magnetic" flux enclosed by the ring. 

decoupled. Here, the propagation of the atoms in the ring 
are symmetric so that anti-symmetric combination does 
not affect the propagation. As a result, we can obtain the 
inhomogeneous TL liquid as the effective model of 
the double Y junction systems, where the TL parameter 
on the left and right lead is g, and that of the symmetric 
combination S is 2g. Following Ref. [TtJ, we can investi- 
gate the transmission and the reflection at the junctions. 
The solution of the equation of motion naturally gives the 
contribution of multiple reflections. Interestingly, the ef- 
fect of normal and Andreev- type reflections cancel out, 
so that the total transmission approaches 1 as the elapsed 
time approaches infinity. Experimentally, real-time ob- 
servations of the density could resolve each reflection. 

Now let us consider the effect of a "magnetic flux" cf> 
inside the ring. The present geometry is a typical setting 
to observe the AB effect. If we consider the transport 
of a single charged particle in the same geometry with 
an empty ring, the transmission probability from the left 
lead to the right vanishes when <f> is the half unit flux 
quantum. This is because of the cancellation between 
the probability amplitudes coming from the two paths. 
Although a magnetic flux actually does not induce the 
AB effect for neutral atoms, several possible approaches 
for introducing gauge fields and AB type phase factors 
for neutral atoms have already been discussed in the lit- 
erature [lij]. A common ingredient of many methods is 
that of the Berry's phase arising from the orbital motion 
of atoms [nl . 

We find that, however, in our problem of filled one- 
dimensional BEC, the transport between the leads is 
completely insensitive to the effective magnetic flux </>, 
in the low-energy limit. Namely, the AB effect is absent. 
This rather surprising conclusion follows quite naturally 
from the present formulation, as follows. 

With a gauge transformation, the effect of the flux 
can be described by the boundary conditions tp u (—L) = 
Lpi(—L) at the left junction and ip u (L) = tpi(L) + cj> at the 
right one. They determine the boundary values of the an- 
tisymmetric combination as ip a {—L) = and (p a {L) = 4>. 
The current on the ring is just a superposition of the 
symmetric flow between the two leads described by ip s 
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and the persistent current described by ip a . Only the 
latter depends on the flux (f>. As it is only the symmetric 
combination ip g that couples with the leads, the flux <p 
has no effect on the transmission between the leads. 

If we consider the system with fermions instead of 
bosons in the same geometry, the AB effect must be 
present, at least when the fermions are non-interacting. 
As the non-interacting fermions in one dimension corre- 
sponds to the TL liquid with g — 1, our result might 
appear in contradiction. However, there is actually no 
contradiction, because our analysis applies only to the 
"strongly coupled" b.c. ([6]). In the system of non- 
interacting electrons, the junction is described by a dif- 
ferent b.c. The "strongly coupled" b.c. is expected to be 
stable even in a fermionic system, when the interaction 
is sufficiently attractive (g > 9) [llj. For such a system, 
the AB effect should be absent as well. 

Finally, we discuss the limitations of and possible cor- 
rections to the present results. First, the TL liquid de- 
scription is only valid in the low-energy limit. This re- 
quires T, 2irv I L <C p, where L is the system size (length 
of the branches), and the chemical potential p = U p^ 
gives the ultraviolet cutoff. [3| When a packet of extra 
particles is injected, the variation Ap of the particle den- 
sity corresponds to a variation of the local chemical po- 
tential Ap = UAp. Thus Ap <C p^ is required for the 
TL liquid description to apply. For a larger amplitude, 
the single-particle physics is expected to show up, instead 
of the collective nature characterized by the TL liquid. 

The "strongly coupled" b.c. ((6]), which we assumed in 
the present Letter, is exact only in the low-energy limit, 
even when the bulk is completely described by the TL 
liquid. There are corrections due to finite energy scales 
present in a realistic system. The leading correction to 
the results based on the b.c. (|6|) can be related to, at a suf- 
ficiently low energy scale, the leading irrelevant operator 
with scaling dimension 4g/3. It would give, for example, 
a positive reflection proportional to e 8s / 3 ~ 2 on top of the 
universal negative reflection —1/3. Here e denotes the 
largest among the possible infrared cutoff scales, such as 
T, 2irv/L, and UAp. When the interaction is weak, it is 
difficult to satisfy the conditions for the TL liquid theory 
to be applicable in the bulk. Interestingly, however, the 
corrections at the junction vanishes quickly because g is 
large in such a case. 
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